ARPES spectral functions and Fermi surface for Lai 86 Sr .i4CuO4 compared with 

LDA+DMFT+E k calculations 
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Slightly underdoped high-T c system Lai.geSro.wCiiC^ (LSCO) is studied by means of high energy 
high resolution angular resolved photoemission spectroscopy (ARPES) and combined computational 
scheme LDA+DMFT+Ek- Corresponding one band Hubbard model is solved via dynamical mean- 
field theory (DMFT), while model parameters needed are obtained from first principles within 
local density approximation (LDA). An "external" k-dependent self-energy Ek describes interaction 
of correlated electrons with antiferromagnetic (AFM) pseudogap fluctuations. Experimental and 
theoretical data clearly show "destruction" of the LSCO Fermi surface in the vicinity of the (tt,0) 
point and formation of "Fermi arcs" in the nodal directions. ARPES energy distribution curves 
(EDC) as well as momentum distribution curves (MDC) demonstrate deviation of the quasiparticle 
band from the Fermi level around (ir,0) point. The same behavior of spectral functions follows from 
theoretical calculations suggesting AFM origin of the pseudogap state. 

PACS numbers: 74.72.-h; 74.20.-z; 74.25. Jb; 31.15.A- 



I. INTRODUCTION 

One of the not yet solved puzzles of cupratc high- 
temperature superconductors (HTSC) is a nature of un- 
derdoped normal state - the pseudogap regime 1 . Per- 
haps most powerful experimental tool to access electronic 
properties of the pseudogap state is angular resolved pho- 
toemission spectroscopy (ARPES) 2,3,4 . Common under- 
standing is the fluctuating origin of the pseudogap state, 
however type of the fluctuations is still under discus- 
sion. Whether these are superconducting fluctuations 5 or 
some order parameter fluctuations (AFM(SDW), CDW, 
stripes, etc.) 6,7 coexisting or competing with Cooper 
pairing is up to now undecided. 

There are several prototype compounds among high- 
Tc systems e.g. hole doped Bi 2 Sr 2 CaCu20 8 _ ( 5 (Bi2212) 
system or the electron doped one - Nd2- a; Ce 2: Cu04 
(NCCO). Now a lot of experimental ARPES data on 
Bi2212 and NCCO is available (see the review Ref. 2). 
For instance Fermi surface (FS) maps, quasiparticle band 
dispersions and even self-energy lincshapes mapped onto 
some models are obtained from modern ARPES measure- 
ments 6 . Into this list of prototype componds should be 
of course included the first ever high-T c hole doped sys- 
tem La2- 2; Sr 2 ,Cu04 (LSCO) which was also investigated 
both in theory and experiment with great details. 2 

Within the normal underdoped phase (pseudogap 
regime) a number of interesting physical phenomena were 
discovered. For example, FS is partially "destructed" in 
the vicinity of the so-called "hot-spots" (points of cross- 
ing between FS and AFM umklapp surface). "Shadow 
bands" (partial folding of band dispersion) appear ap- 
parently as a result of short range AFM order. Forma- 
tion of the so called Fermi "arcs" around Brillouin zone 
(BZ) diagonal reminiscent of the parts of noninteract- 
ing FS is experimentally detected in numerous ARPES 



experiments 2 . Despite apparently the same physics be- 
hind, pseudogap regime demonstrates some material spe- 
cific features. For Bi2212 Fermi "arcs" go almost up to 
BZ border where they are strongly blurred. The NCCO 
also has Fermi "arcs" but FS "destruction" looks differ- 
ent. The "hot-spots" are well observed in NCCO while 
towards the BZ border FS is almost restored as a non 
interacting one. 8 

The present paper is devoted to pseudogap behavior in 
the underdoped LSCO and its comparison with Bi2212 
and NCCO. 

According to common knowledge high-T c systems 
are usually doped Mott insulators, effectively described 
by the Hubbard model. Most common method to 
solve the Hubbard model in our days is dynamical 
mean-field theory (DMFT) 9 . Its exactness in the infi- 
nite spatial dimensions limit makes DMFT a local ap- 
proach. Since high-T c compounds as it is well estab- 
lished have quasi two-dimensional nature spatial fluctua- 
tions play important role for their physics. To overcome 
this difficulty we introduced DMFT+Sk computational 
scheme 10,11,12 which supplies conventional DMFT with 
"external" k-dependent self-energy. Main assumption of 
the DMFT+Sk is additive form of the self-energy which 
allows one to keep conventional DMFT self-consistent 
set of equations. The DMFT+Sk approach was used 
to address the pseudogap problem 11 , electron-phonon 
coupling in strongly correlated systems 13 and disorder 
induced metal-insulator transition within the Hubbard- 
Anderson model 14 . For the pseudogap state this self- 
energy Sk describes the interaction of correlated elec- 
trons with non-local (quasi) static short-ranged collec- 
tive Heisenberg-like antiferromagnetic (AFM or SDW- 
like) spin fluctuations 15 . DMFT+Sk approximation was 
also shown to be appropriate to describe two-particle 
properties e.g. optical conductivity. 16 
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As a possible way of theoretical simulation of 
the pscudogap regime for real materials we proposed 
LDA+DMFT+E k hybrid method 7 . It combines first 
principle one-electron density functional theory calcula- 
tions within local density approximation (DFT/LDA) 17 
with DMFT+E k . 18 

LDA+DMFT+E k method allowed us to obtain Fermi 
arcs and "hot-spots" behavior for both electron doped 
e.g. Ndi 85 Cc i 5 Cu0 4 8 (NCCO) and Pr x 85 Ce i 5 Cu0 4 
(PCCO) 19 as well as hole doped Bi 2 Sr 2 CaCu 2 8 _« 
(Bi2212) 7 high-T c cuprates. Pseudogap behavior of dy- 
namic optical conductivity within LDA+DMFT+Ek 16 
was also discussed for B12212 7 and NCCO 8 . 

This communication reports LDA+DMFT+Ek com- 
putations of Fermi surface and spectral functions for hole 
undcrdoped Lai.sgSro.MCuC^ (LSCO) system supported 
by high energy, high resolution bulk sensitive ARPES 3 . 

II. COMPUTATIONAL DETAILS 

The La 2 Cu04 system has base-centered orthorhom- 
bic crystal structure with space group Bmab with two 
formula units per cell 20 . Corresponding lattice parame- 
ters are a=5.3346, a=5.4148 and c=13.1172 A. Atomic 
positions are the following: La(0. 0,-0. 0083, 0.3616), 
Cu(0,0,0), 0(1/4,1/4,-0.0084), 02(0.0,0.0404,0.1837). 

As a first step of LDA+DMFT+Ek method we per- 
formed density functional theory calculations in the local 
density approximation (LDA) for these crystallographic 
data. Band structure was obtained within the linearized 
muffin-tin orbitals (LMTO) method 21 . It is well known 
for these compounds that Fermi level is crossed by an- 
tibonding 02p-Cu3d partially filled orbital of x 2 — y 2 
symmetry. Tight-binding parameters for this band were 
calculated by N-th order LMTO (NMTO) method 22 as 
(in cV units) t = -0.476, t' = 0.077, t" = -0.025, 
t'" = —0.015. These values agree well with previous 
studies. 23 Coulomb interaction value on effective Cu- 
3d(x 2 — y 2 ) orbital was calculated by constrained LDA 
approach 24 and was found to be £7=1.1 eV. These LDA 
obtained parameters are used to set up corresponding 
one-band Hubbard model. 

The second step is further treatment of the above de- 
fined Hubbard model within the dynamical mean-field 
theory (DMFT) self-consistent set of equations 9 supplied 
by "external" momentum-dependent self-energy Ek 11 . 
Using the additive form of self-energy (main approxima- 
tion of the scheme neglecting the interference between 
Hubbard interaction and pscudogap fluctuations, which 
allows one to preserve conventional DMFT equations) 
one can define LDA+DMFT+Ek Green function as: 

Gk ^ = ~L m vTl (1) 

oj + fj, — e(k) — E(cj) — Ek(o;) 

where bare electron dispersion e(k) is defined by LDA 
calculated hopping parameters listed above. To calcu- 
late Ek we used a two-dimensional pscudogap model 1,15 



describing nonlocal correlations induced by (quasi) static 
short-range collective Heisenbcrg-likc AFM spin fluctua- 
tions. Thus we introduce correlation length dependence 
of the pseudogap fluctuations into conventional DMFT 
loop. 

There are two points which make DMFT+Ek different 
from the usual DMFT. First, momentum dependent Ek 
is recalculated on each DMFT iteration (Ek(/i, oj, [E(w)]) 
is in fact the function of DMFT chemical potential and 
DMFT self-energy). Second, DMFT+E k lattice problem 
is defined at each DMFT iteration as: 

= N ^ u + n - e(k) - S(w) - EkH ' (2) 

After numerical self-consistency is reached we get the 
Green's function (1) with corresponding E(w) and E k (w) 
taken on the last DMFT iteration. All further computa- 
tional details can be found e.g. in Refs. 7,8,11. 

As an "impurity solver" for DMFT equations numer- 
ical renormalization group (NRG 25,26 ) was employed. 
Temperature of DMFT(NRG) computations was taken to 
be 0.011 eV and electron concentration used was n=0.86. 

Self-energy Ek(w) due to pseudogap fluctuations de- 
pends, in general, on two parameters: the pseudogap 
amplitude A, and the correlation length £. 1,15 The value 
of A was calculated as in Rcf. 11 

A 2 = u 2 <n^n il > K ( )2>; (3) 

where local densities n^, and double occupancy 
< njtriji > were calculated within the standard DMFT 9 . 
Behavior of A MS M function of hopping integrals and 
Coulomb interaction was studied in our previous work 11 , 
while A as a function of occupancy n was investigated in 
Ref. 7. For £ we believe it is more safe to take experimen- 
tal values. In this work the value of A was calculated to 
be 0.275 eV and £ was taken to be 10a, where a - lattice 
constant. 27 



III. EXPERIMENTAL DETAILS 

The high-energy ARPES measurements were carried 
out at BL25SU in SPring- 8, using incident photons of 
500 eV, on single crystal samples. The normal to the 
cleaved sample surface was set almost parallel to the axis 
of the lens of the analyzer and the sample was set to 
about 45° from the incident light direction. The pho- 
toelectrons within polar angles of about ±6° around the 
normal to the sample were simultaneously collected using 
a GAMMADATASCIENTA SES200 analyzer, thereby 
covering more than a whole Brillouin zone along the di- 
rections of the analyzer slit. The Fermi surface mapping 
was performed by changing the angle along the perpen- 
dicular direction to the analyzer slit. The base pres- 
sure was about 4 x 10 _8 Pa. The (001) clean surface 
was obtained by cleaving the samples in situ in vacuum 



at the measuring temperature of 20 K. The overall en- 
ergy resolution was set to 100 and 170 meV for high- 
resolution measurements and Fermi surface mapping, re- 
spectively. The angular resolution was ±0.1 (±0.15) for 
the perpendicular (parallel) direction to the analyzer slit. 
These values correspond to the momentum resolution of 
±0.0247r/a ( ±0.0367r/a) at W=500 eV, where a is twice 
the Cu-0 bondlcngth within the CuC>2 plane. By virtue 
of the longer photoelectron mean free path of ~12 A at 
the kinetic energy of ~500 eV than that for conventional 
ARPES at hv -20-60 eV, the bulk contribution to the 
spectral weight is estimated as about 60%. The position 
of the Fermi level was calibrated with Pd spectra. 



IV. RESULTS AND DISCUSSION 
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FIG. 1: LCOO Cu-3d(x — y ) band dispersion along high 
symmetry directions of square Brillouin zone computed with 
LDA+DMFT+E fe . The Fermi level is zero. 

In case of finite temperature and interaction strength 
we have to take into account the finite life time effects 
of quasiparticles. Thus instead of just a dispersions e(k) 
one should rather work with the spectral function A{lj 1 k) 
given by 



A(u,k) = --ImG(cj,k), 

7T 



(4) 



with retarded Green's function £?(w, k) obtained via 
LDA+DMFT+Sk scheme 10 - 11 ' 16 . Of course there are 
considerable lifetime effects originating from the £k cor- 
responding to interaction with AFM fluctuations (substi- 
tuted in our approach by the quenched random field). 

In Fig. 1 contour map of spectral function (4) obtained 
from LDA+DMFT+S k for Cu-3d(x 2 - y 2 ) band is pre- 
sented. The width of the spectral function is inversely 
proportional to the lifetime. Around (7r,0) point one 
can clearly see splitting of the spectra by AFM pseu- 
dogap fluctuations of the order of 2 A. Also AFM nature 
of the pseudogap fluctuations leads to formation of the 
"shadow" band which is much weaker in intensity and 
becoming the real quasiparticle band in case of complete 
folding in case of long-range AFM order. 

Figure 2 displays experimental energy distribution 
curves (EDC) on the panel (a) along (0,0)-(7r,0) direc- 




:b): 
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FIG. 2: ARPES EDC curves (a) and LDA+DMFT+E fe spec- 
tral functions (b) along (0,0)-(7T, 0) high symmetry direction; 
(c) ARPES MDC curves around (tt, 0) point for LSCO at 
x=0.14. On panels (a), (b) and (c) filled circles guide the 
motion of the A(uj,\C) maxima. The Fermi level is zero. 



tion. Around (7r,0)-point as shown by stars certain de- 
viation of the A(w,k) maxima from the Fermi level a 
kind of "turn-back" is observed. We attribute such be- 
havior of the A(w,k) to pseudogap fluctuations. Similar 
theoretical behavior is shown on panel (b) 28 of Fig. 2 as 
calculated by our LDA±DMFT±£k approach (see also 
Fig. 1). The same behavior is also observed as traced by 
circles in experimental ARPES momentum distribution 
curves (MDC) demonstrated on panel (c) of Fig. 2. 

The bulk-sensitive high-energy ARPES data for 
Lai.86Srn.i4Cu04 show a clear "turn-back" structure of 
the EDC peak as a function of momentum near (0,— 7r), 
which were not seen in the previous low-energy ARPES 
data for Lai.ssSrn.isCuO^ 4 . The contour map of the 
spectral weight in the vicinity of Ep seems to be essen- 
tially similar in overall features for this doping level be- 
tween the high-energy and low-energy ARPES studies. 

Experimental and theoretical Fermi surface maps are 
shown in Fig. 3 at panels (a) and (b) correspondingly. 
Both pictures reveal strong scattering around (7r,0)-point 
which we associate with scattering in the vicinity of the 
so-called "hot-spots" (cross-points of a Fermi and AFM 
umklapp surfaces) which are close to the (7r,0) 7,8 . Such 
strong scattering comes from scattering processes with 
momentum transfer of the order of Q=(n, 7r) 1,15 , corre- 
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FIG. 3: Fermi surfaces of LSCO at x=0.14 from experiment 
(left panel) and LDA+DMFT+Ek computations (right panel) 
Red crosses on the left panel correspond to experimental ki? 
values. 
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sponding to AFM pseudogap fluctuations. Along nodal 
directions we observe typical Fermi arcs. They are pretty 
well seen in the theoretical data while for experiment we 
have just narrow traces of them. 

V. CONCLUSION 

Our LDA+DMFT+Sk hybrid approach was shown to 
be an effective numerical tool to describe short-range or- 
dered state in quasi-two dimensional systems. 7,8,19 Mate- 
rial specific model parameters such as hopping integrals 
(which define bare electronic band dispersion of effective 
Cu3d(a; 2 — y 2 ) orbital) were calculated via LDA based 
NMTO method. 22 Coulomb interaction parameter U was 
obtained from constrained LDA method. Pseudogap am- 
plitude A was computed using LDA+DMFT scheme. 8,11 
Supplementing conventional DMFT self-energy by Sk(^) 
describes non-local dynamical correlations due to short 
ranged collective Heisenberg like AFM spin fluctuations. 

In this work we performed LDA+DMFT+S k calcu- 
lations for hole-doped Lai.g4Sro.i4Cu04 compound in 
the pseudogap regime. Because of fluctuations of AFM 
short-range order we clearly observe formation of the so 
called "shadow bands" as partially folded bare dispersion. 
Pseudogap is formed around (tt,0) point which is quali- 
tatively the same as in Bi2212 7 , NCCO 8 and PCCO 19 . 
As to Fermi surface of LSCO it is alike that obtained for 
Bi22I2 7 . Namely the "hot spots" are not well resolved 



since the crossing point of the bare Fermi surface and 
AFM umklapp surface are very close to Brillouin zone 
border. This fact is essentially due to the shape and size 
of LDA Fermi surface. In this respect, situation here is 
different from that for NCCO 8 and PCCO 19 , where "hot 
spots" are clearly seen. To support these theoretical re- 
sults we present here the new high energy, high resolution 
ARPES data for LSCO. Typical pseudogap- like effects of 
Fermi surface destruction were observed in both theory 
and experiment. The same is true for spectral functions. 
The overall semiquantitative agreement between theory 
and experiments basically supports our general picture of 
the pseudogap state as due to strong scattering of carriers 
by short-range AFM order fluctuations. 
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